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1 Introduction and main result 



CM. 

' This paper is devoted to study the stabihzation of the Unear wave equation in a bounded domain 

(— I , damped in a subdomain when the geometrical control condition (see | ELRp of the work of C. Bardos, 

G. Lebeau and J. Ranch is not fulfilled. In such case, they | BLH,| proved that the uniform exponential 
decay rate of the energy cannot be hoped due to the existence of a trapped ray that never reaches the 
support of the damping. Another important contribution in this field was done by G. Lebeau | Lej who 
establishes a logarithmic decay rate for the dissipative wave equation when no assumption on rays of 
04 ' geometrical optics is required, but when more regularity on the initial data is allowed. Now, it seems 

<*' , natural to search a general description of the geometries of both domain and support of the damping 

under which the energy of the dissipative wave equation decreases in a polynomial way. A first answer 
in this direction was done by Z. Liu and B. Rao \ LR) who consider the wave equation on a square 
' damped in a vertical strip. In this paper, we improve the geometry to a partially cubic domain where 

the damping acts in a neighborhood of the boundary except between a pair of parallel square face of 
CO ■ the cube. 

o 

fH , Before stating the main result of this paper, we begin by presenting precisely the geometry of our 

' problem. Next, we introduce the equations that will be used throughout this work. Our main result 

g . is given at the end of this section. 

> 

I 1.1 The geometry 

H , 

^ Let mi, TO2, p > 0. Let f2 be a connected domain in R'^ bounded by Fi, T where 

Fi = [—mi, nil] X [—1712, x with boundary dVi, 
F2 — [—mi, mi] X [—m2,m2\ x {—p}, with boundary 9F2, 

T is a surface such that T C \((— rrii, mi) x (— m2, m2)) x R, with boundary dT — dVi U 9F2. 

We suppose that the boundary of il, dfl = Fi U F2 U T, is C°°. Let be a neighborhood of T in R'^ 
and let oj — finQ. 

As is a neighborhood of T in M.'^, there exists E (0, min (mi, m2, p) /2) such that [x^i — 2ro, Xbi + 2ro] x 
[xb2 - 2ro, Xb2 + 2ro]x[xb3 ~ 2ro,Xbz + 2ro] C 8 for any {xbi,Xb2-,Xb3) G dTi. Let ho = min ^1, (ro/8)^^ . 
Next, we choose — (— + ro, mi — To) x (— m2 + ro,m2 — To) x (— f , |)- 

Throughout this paper, c denotes a positive constant which only may depend on (mi, m2, p). Also 
7 will denote an absolute constant larger than one. The value of c > and 7 > 1 may change from 
line to line. 
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1.2 The equations 



We consider the dissipative wave equation in O with initial data {wq, wi) G H"^ (O) fl {Q) x Hq {Q). 

dfW - Aw + a (.r) diw =0 in O x K+ 

w ^0 ondnxR+ (1.1) 

{w{-,0) ,dtw{-,0)) = {wo,wi) in , 

with a non-negative dissipative potential a G L°° (O) such that a > in w. Denote 

£{w,t)= f (\dtw{x,t)\^ + \Vw{x,t)\^)dx = S{w,0)-2 [ [ a{x)\dtw {x,0)\^ dxd0 , 
Jq^ ' Jo Jn 

and recall that £ {w, t) is a continuous decreasing function of time. 

Denote by {lJ'j}j>i, < /xi < /i2 < Ats < • • •, the eigenvalues of —A on O with Dirichlet boundary 
conditions and by {f.j}j>i the corresponding normalized eigenfunctions, i.e., ||^j||i,2(Q) = 1. Let u = 
u{x,t) be the solution of the wave equation 

d^u - Au =0 in X R 

u =0 ondnxR (1.2) 

{u{-,0) ,dtu{-,0)) ={uo,ui) in O . 

Suppose that uo= Yl ^^'^ = X] bjlj are such that 

l|wo||^2nHi(f^) =J2f^^3 l^?r < and ||ui||Hi(a) =Y,f^o I^^T < +0° . 
then it is known that 

and u e C (M; (fi) n iji (O)) n (M; (Q)) n (R; (fi)) . Let 

g{u,t) = j {\dtu{x,t)\^ + \Wu{x,t)\^^dx = g{u,Q) . 



b" cos (ty/JIj) + 



■ sin (i. 



1.3 The bicharacteristics 

It is usual to associate with the wave equation, the geodesies. Recall that the bicharacteristics are 
curves in the space-time variables and their Fourier variables described by 

x(s) =Xo + 2^{s)s / ^(s) =^o 



t (s) = to — 2t {s) S \ T (s) = T( 

with \^{s)f - r2(s) = for s e [0,+oo), when (xo, ^o, r^) e IR''xM^\{0}. The rays are the 
projection of the bicharacteristics on the space-time domain 

x (s) — Xo — 2^oS = 
t (s) + 2toS = 
\^of-r^ =0, 

here, to = and Tq ^ 0. The generalized rays are rays taking into account the geometry by following 
the rules of optic geometric. 
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The key geometric observation in our setting is that we do not have that any generahzed ray meets 
Lo. More precisely, any [xo^ to, Co, To) G i^o x R x \{0} such that = (0, 0, ±1) generates a trapped 
generahzed ray which never goes outside [—mi, mi] x [— m2, TO2] x [~P, p\- As a result, we do not have 
an uniform exponential decay for the dissipative wave equation for any damping only acting in w. Of 
course, the logarithmic decay rate still holds from Carleman inequalities. However, we may hope a 
better decay rate because our trapped generalized ray behaves quite simply by bouncing between Fi 
and r2 always in the same direction = (0, 0, ±1). Furthermore, observe that the geometrical control 
condition is fulfilled with w U tt>o . 



1.4 Main result 



Now we are ready to formulate our main result. 

Theorem.- There exist C > and 5 > Q such that for any t > and any initial data (wq, wi) G 
H'^ (il) n Hi (il) X Hi (fi), the solution w 0/ (IF1() satisfies, 

(|(9tw(a;,i)|V |Vu)(x,t)|^) da; < ^ ||(wo,wi)||^2nffi(n)xffi(n) • (1-3) 



Our strategy to get such polynomial decay rate will consist to establish an observability estimate 
for the wave equation u solution of (|1.2|) . or more precisely, an inequality which traduces the unique 
continuation property for u between u x (0, T) and x {0} for some T < +00. For example, it is now 
known that an interpolation inequality of the form 

, t) I dxdt 



where A = — rr-. , imphes 11.311 . On another hand, it is not difficult to deduce from 

Il(«0'"i)llffi(n)xi2(n) 

an inequality in the following form 

ll("o,ui)ll_f/i(f2)xL2(o) ^ j^a{x)\dtu{x,t)\ dxdt 

and conversely, the later inequality implies where the constant C may change of values. In this 
paper, H1.3|l comes from 

/•(CA)" , 

|2 ^ ^ AT / / 10 / iM2 



("o,Mi)|lHi(n)xL2(ii) < C' / / \dtu{x,t)\ dxdt 



As it was suggested in the previous subsection, we have to pay more attention on a ray of geometrical 
optics bouncing up and down infinitely between the two parallel planes Fi and F2. In same time, we 
need to estimate in a good way the dissipation phenomena in order to improve the logarithmic decay 
rate. To this end, we apply some simple tools usually used in the propagation of singularities | AG| - 
I CV| -|~n | in order to link the number of reflections and the s variable of the bicharacteristic flow. In 
particular, we will work with the operator Op (D) — ids + ft. (A — df) for h G (0, ho]- Observe that 
the product of four, mono-dimensional, solutions of the Schrodinger equation ids i hd^ can create a 
solution of Op (D) a {x, t, s) = 0. The dispersive property of the linear Schrodinger equation ids ± hd^ 
will be exploited. 

The next section describes an interpolation inequality. In section 3, we give the proof of Theorem. 
We close this paper with two Appendixes devoted to prove some technical results. 



3 



2 Interpolation inequality 



The purpose of this section is to estabhsh the following inequality. 

Theorem 2 .- Let T > 0. There exist C > and 7 > 1 such that for any h e (0, ho] and initial 
data (uo, ui) G H'^ (fi) D Hq (fi) x Hq (Q), the solution u of p.2|) satisfies 

1 /2 

+C\//l ||(wo,Ui)||^2n//i(o)xHi(a) ll("0,Ul)||_ffi(f2)xL2(O) • 

The rest of this section is devoted to the proof of Theorem 2. We begin to introduce some weight 
functions inspired by with the properties of localization, propagation and dispersion. Next, we 
make appear the Fourier variables and introduce some Fourier integral operators dependent on the 
number of reflections between Fi and F2. Their properties at the boundary are analyzed. The process 
of propagation is then applied. Some parameters are adequately chosen and that will complete the 
proof. 



2.1 The weight functions 

The hypothesis saying that T C \((— mi, mi) x (— m2,m2)) x M implies that for any Xq G uJ^, 
B {Xo,ro/2) n d^l = 0, where B{xo,r) denotes the ball of center Xg and radius r. We introduce 
Xxa = X G Cq" (^o, ^0/2)) be such that < x < 1 and x = 1 on B {xo, fo/4). 

From any point Xq G oJ^, we will localize around Xo and eventually be able to propagate some local 
regularity. To this end, let /i G (0, ho] and let us define 

X / 1 __L^\ / 1 1 

a {x, t,s) ^ \ — e 4h ..+1 e -"-=+1 

Oo (x, t) = a{x — Xo, t, 0) and ip [x, t) = x {'^) '^o {x, t) . 
We get the following identities 

{ids + h{A~ dj)) a {x, t,s) = V {x,t,s) e Q x R x (0, +00) , 
1 1 



\aix,t,s)\^ e . ji^^ _ (2.i) 



Now we use such weight functions Oq and (p with u the solution of the wave equation p.2|l as follows. 
By integrations by parts, 

InxRXix) \aodtu{x,t)f dxdt 
= ~/nxR^(^)^* (l"o (^'*)l^) 5^* (l"(2;>OI^) dxdt - J^^^ x (x) \ao {x,t)\'^ dfu {x,t) u {x,t) dxdt 
= /f2xR^(^) h^t (ko ix,t)\'^j \u{x,t)\^ dxdt - /oxRX(a;) 1^0 ix,t)fdfu{x,t)u{x,t)dxdt . 



As d; 



(^\ao {x, = — jflo {x, + t"^ jflo (x, i)|^ and t'^ \ao {x, t)|^ < 4 [oo (a;, i/v^ |^, we have 

/nxR^(^) |ao9fW (x,t)|^ dxdi 
- ^ /oxk'^o (a^j*/V2) (/3 (a;,i/V2) \u {x,t)\'^ dxdt + /q^rOo {x,t) (p {x,t) dfu {x,t) u {x,t) dxdt 



(2.2) 
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Let {Gi; i € /} be a family of open sets covering uJq, i-e., C IJ Gi, such that Gi — <\x — x 

where {2:0}^^^- € TJ^ and / is a countable set such that the number of elements of / is bounded by -j^^ 
for some constant Cq > independent of /i £ (0, /iq]. Consequently, 

L„x(o,T) \dtu{x,t)\^ dxdt < e5^'X.„x(o,T)S~^*' \dtu{x,t)\^ dxdt 

< +2 ^/^^^jjXi:^ (2;) |a(a;-2:o,i,0)atu(x,i)| dxdt ^23) 

< e5^'+2 X^/nxR^a;' (2;) \a{x- xl,t,{))dtu{x,t)\^ dxdt 

and we will search to bound /qxrXx» (2;) \a (x — xl,t,(fj dtu {x,t)\^ dxdt by a suitable term Eh{u) 
independent of i in order to get 

\dtu{x,t)\'dxdt< ^°-^^(") 
WoX(o,T) /ly/i 

for some Co > independent of u and h. To this end, we will first study in the next subsections the 
second term of the second member of (|2.2|) , 

Co (x, t) (f [x, t) f {x, t) u (x, t) dxdt when / = ^^u. 



2.2 The Fourier variables 

Denote 

^ , t)= [ e-^(^«+*^) ^ (x, t) f [x, t) dxdt , 
JnxE 

then for any (x, i) e x M, 

flo (x, i) ^ (x, t) / (x, t) = ao (x, t) / e^(-«+*^) ^/ (C, r) d^rfr . 

(27rj JR4 

Let A > 1. We cut the integral over r e M into two parts, {|r| > A} and {|t| < A}. Next, for {|t| < A} 
and ^ = (^1, ^2,^3)1 the integral over ^3 e R is divided as follows. Denote (2Z + 1) = {2n + 1 \ri e Z}, 
then 

ao (x, t) ip (x, t) f (x, t) = a, (x, t) ^ E /r. //^-i' /|.|<a (C, ^^^r 

+«o (X, t) ^ /j,3 /|.|>, e^(^«+*^' (C, r) didr . 

Consequently, 

/nxR ^) ^) " ^) '^^'^^ ~ ^0 

= /oxM «o (5^)- E /r^ /|.|<A ?! it r) d^dr u (x, t) dxdt (2-4) 

«o3e(2Z+i) 

where 

i?o = / ao (x, i) f [ e^(-«+*^) ^ r) d^dr u (x, t) dxdt . 

JnxM (27r) Jr3J|t|>a 

Moreover, it holds 

-Ro < c ^J^y/G {u,0)VG {dtu,0) . (2.5) 
The proof of ()2.5|1 is given in Appendix A. 
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2.3 The Fourier integral operators 



Denote x = {xi,X2,X3) e Jl, f = (6,676) with (6,6) ^ K^. Let (^0,63) = {xoi, Xo2, Xo3, ^03) G 
uJ^ X (2Z + 1), we introduce for all s > and n e Z 

(x„,63)/(a;,i) 

a f .Ti — Xoi — 26^5, — Xo2 — 26^5, 2^3 " ("1)" + + 26/is , ^ + 2t/is, s j dr 



Let (P, Q) e N^. Consider the solution 



2P+1 



^,p,Q (xo, 63) / (a^, t) = 51 ^?(a^o,63)/(a;,i) . 

One can check that for any {xo, 63, P,Q) x (2Z + 1) x N^, 

{ids + h{A- d^)) As^p^Q {xo, 63) / (x, t) = V (x, s) e X R X (0, +c3o) 



(2.6) 



2.4 At s = 



Since 



2^ r /-Sos+i 
^0 (2^0,63)/ (a;,i) Oo (x,i) — — 4 / / 



e^(.?+tr) ^(^^^) ^^^^^ 



we then obtain from (|2.4|) that 

o {x,t) (p {x,t) f {x,t) u {x,t) dxdt — Ro = / >^ (xo, 63) / (2^, ^) u{x,t)dxdt 



On 

/nxR 
Observe that 



?o3e(2Z+i) 



2P+1 



^o,p,g (xo, 63) / (a;, = ^0 (a^o, 63) / {x, t)+J2 ^0 i^o, 63) / {x, t)+ ^ {xo, 63) / (a:, t) , 

n=l -20<n,<-l 

with the convention that for Q = 0, ^ Aq {xq, 63) / (a^, t) — 0- So, we deduce that 

-2Q<ri<-l 

/ aoix,t)if{x,t) f {x,t)u{x,t)dxdt-Ra-Ri = ^ A^^p^q (xo^^oa) f {x,t) u {x,t) dxdt 

(2.7) 



4o36(2Z+l) 



where 

-Ri = - 



OxR 



E 



:o3e(2Z+i) 



2P+1 



E ^0 (^o, 63) / {X, t)+ ^0 (^o, 63) / {X, t) 

n=l -2Q<?1<-1 



u(x,t) dxdt 



We estimate i?i uniformly with respect to (P, Q) as follows. 

^1 </axR E E 1^^(2:0,63) /(a;, 01 \u{x,t)\dxdt 
C„3G(2Z+i) Tiez\{o} 

</nxR E (Sif /r^ i|r|<A "^Z ^) '^^^^ 

C„3e(2Z+i) ^ ^ 

X) aixi- Xoi,X2- Xo2,{-'^T X3 + 2nA^p- Xo3,t,Q] \u{x,t)\dxdt 
nez\{o} ^ ^ 

^ (2i)^/R3i|r|<A|^/(^'^)| '^^^^ 

/nxR a{xi-Xoi,X2-Xo2-X~'^T X3 + 2np-Xo3,t,Q) \u{x,t)\dxdt 

nez\{o} 
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a (^Xi - Xol,X2 - Xo2, {-ir X3 + 2n||^p - Xo3,t,(ij 

~Th (^i-^oi)^ + (^2-a;o2)^ + ((-l)"2:3+2ri^s3^p-a;o3) 



Now, notice that 



e * . 

The hypothesis saying that T C M^\((— mi,mi) x (— m2,m2)) x R imphes that for any Xo G uJ^ 
and X — {xi, 0:2, 2:3) S Q, such that {xi, X2) ^ [—mi + ro/2, mi ~ ro/2] x [— TO2 + ro/2, m2 — '"o/2], we 
have (xi — a;oi)^ + (2^2 ~ 2:02)^ ^ (''o/2)^, but such hypothesis also imphes that for any Xg G uJ^ and 
X = (xi, X2, X3) e such that (xi, X2) G [— m-i + ro/2, mi — ro/2] x [— m2 + ro/2, TO2 — ro/2], we get 
X3 e [— p, p] and therefore (—1)" 2:3 + 2n-^^p — Xg^ > jp for any n lE Z \{0} . So, for any n e Z \{0} 



a \^xi - Xoi,X2 - Xo2, (-1)"2:3 + 27i-j|2Lp_ a;o3,i,0 ) < e ''e 
It follows that 



-fg-8t((-l)"-3 + 2«^p-.„3)' 4 



/ ^ a (2:1 — Xoi, 2:2 — 2;o2, (— 1)" 2:3 + 2np — 2:03, 0) \u{x,t)\dxdt < c e y/G (u, 0) 

2nXR ,„-, 



■ neZ\{0} 

Now it remains to compute J^^ i|T|<A '''■^ ''') 



|r|<A 



d^dr. We can check that there exists 7 > 1 such that 



The proof of (|2.8|l is given in Appendix A. 
We conclude that 

Ri<c 



(2.9) 



2.5 On the boundary {xs = ±p} 



Since a (2:1, 2:2, 2:3, t, s) = a (xi, X2, —X3, t, s), we get the following identity 

)"t7^P,H Vn e Z . 

(2.10) 



(a^o, Cos) / ( XI, 2:2, (-1)" g^p, < j = -A^+' (x,, f (^2:1, X2, (-1)" g| p, i j Vn e 



Thus, 



{xo,£.o3)f (^xi,X2, il^Lp^t^ ^ (2;o,Co3) / ^2:1,2:2, il^^'^ ) e Z 



so that 

p 



A,,P,Q{Xo,U)f(xi,X2,^^P,t) = E [A2"(2o,Co3)/(2:,t)+A2»+l(2o,Co3)/(2:,t)] 

= Vs > . 

Also, l|2.1U|) implies 

Al"+' {xo,U) f (xi,X2,-^^P,t^ = -Al^+' {xo,U) f (^xi,X2,-^^P,t^ VrieZ, 
therefore 

As,P,Q (2:0,^03)/ (2^1'2^2,-||ff|P,i) 

= A72Q(a,„,^„3)/(a;i,22,-|f^p,t) +Ar+M2^o,eo3)/(2:i,2;2,-^p,i) Vs > . 
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2.6 The key identity 



By multiplying H2.6|l by u{x,t) and integrating by parts over x M x [0,L], we have that for all 
(xo, eZJ^x {2'Z+l), for all (P, Q) £ and all L > 0, 

/nxE ^o,p,Q {xo, £.03) f {x, t) u {x, t) dxdt = J^^^ ^l,rq {xo, £.03) f {x, t) u {x, t) dxdt 

+ih /gQ^jg ^s,p,Q [xo, £03) f (x, t) duu {x, t) dxdtds . 

Consequently, combining the later equality with (|2.7|) . we have the following key identity 

/oxR ^) ^ / ^) " ^) '^^'^^ ^ Ro~ Ri 

C„3e(2Z+i) (2.11) 

+*^/o /anxR A,p{i^s),Q{i^s) (xo, £03) f {x,t)d^u{x,t) dxdtds 

4o3e(2Z+i) 

for any L > and (P (■fos) , Q (C03)) G N^. L will be taken large enough in order that the first term 
in the second member of H2.11|l is polynomially small like and this uniformly with respect to 

i£o3,P {£03) , Q i£o3)) by using the dispersion of lUJ. Next, for each G (2Z + 1), (P (^03) , Q {^03)) 
will be chosen dependent on (^03 , L) in order that the second term in the second member of (|2.11|l is 
exponentially small with respect to h on the boundary Fi U 



2.7 The internal term 



In this subsection, we study the internal term appearing in H2.11|l 

/ ^ ^L,p,Q {xo, £03) f [x, t) u (x, t) dxdt . 

?o3e(2Z+i) 

First, we have a uniform bound with respect to (P, Q) 

InxR E Al.p,q{xo,£,o3) f {x,t)u{x,t)dxdt 



?o3e(2Z+i) 



5o3e(2Z+i) 



2P+1 

E Al{Xo,£o3)fix,t) 
n=-2Q 



u{x,t) dxdt 



(2.12) 



< E E InxR'^LiXo,£o3)f{x,t)u{x,t)dxdt 
4„3e(2Z+l) nGZ 

Recall that 

Al ixo,U).f{x,t) 

g 4/1 + l g 4/1 + 1 g 



1 (t + 2T/lI.)^ 

g 4 -i/ii + l 



I, (x) 



J>1 

^ E (I + 247) (-) + E (I - 4^) e-^*v^ (x) 
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Therefore, we write 



/nxR i^o, / {x, t) u {x, t) dxdt 



1 (xi-a;oi-2{ihI,) 



f e 

On another hand, 



itr 



iL+l 
1 



1 (x2-x„2-2i:2''-L)^ 

e •-f'+i e 



J ((-l)"^3+2n 



«o3 



(2.13) 



^/-^hL^l 



j. (t + 2-rfi,.L)^ 
g 4 -i/ii + l 



dt dr 



4 -ihi+i dt 



where we have used the foUowing formula 



(2.14) 



e 2 



(r) 



Re z > and z 7^ 



Consequently, by 



LxR (^o' ^«3) / {x, t) u (x, t) dxdt 



(iI+Tp72T2^ 



^?o3-l J|r|<A 

1 (xi-x„i-2;iM,)2 



1 (x2-x„2-2i;2''-'^) 



^ ((-l)"x3 + 2„^p-x„3-2;3>.I.) 



dxdS, 



Thus, 



< 



InxR {xo-, 6) / (a^, *) u {x, t) dxdt 

2V¥ 

e 



dr 



(-l)".3+2n^p-.„3-2C3/iLj ^ 



( 



didT 



Then, 



IfixM^li^oAos) f {x,t)u{x,t)dxdt 
1 



< c 



by using 



(yr^+T) 



3/2 



4 + cv//i(L2 + l) ^ |6" 



i>i 



d£,dT 



(2.15) 



^ e-^(<-^^"-- W-^^3/.i.4^) < 4 + cV /^ (L^ + 1) 
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and \ij {x) \ dx < c ||^j||^2(Q-) — c. On another hand, using Cauchy- Schwartz inequality, we have 

(2.16) 



< 



3>1 



<c^g{dtu,0) 



because 



1 ^ 

E - <^E 



1 



2/3 



< +00 



We conclude from itmil . ^H^-^J^ and lEiHIl, that for any (P, Q) e N^, 



< c 



Al.p,q{xo,S.o3) f {x,t)u{x,t)dxdt 

?„3e(2Z+i) 

^-^) /r3 /|r| 



1+L\/L 
1 /A\T' 



r|<A 



(2.17) 



2.8 The boundary term 



In this subsection, we study the boundary term appearing in (|2.11|l 



ih 



nE ^^,-P(Co3),Q(«o3) i^o, S.03) f [x, t) d„u {x, t) dxdtds 



■«o3e(2Z+i) 

Recall that dVl = Fi U r2 U T. We begin to estimate 

(•L 



ih 



J(riur2)xE 



E ^s.-P(«o3),Q(«o3) (2^0' $03) / (a;, t) d^u (a;, t) dxdtds . 



«o3e(2Z+i) 



First, it holds 

r-L 



lo /(riur,)xR E ^s,p(c<,3).Q(?o3) i^o, ^os) / (a;, t) d^u {x, t) dxdtds 
«o3e(2Z+i) 



Jo J— mi J —mo JM. ^-^ 



«o3e(2Z+i) 



S,P(?o3),Q(«o3) 



Jo J-mi J-m2 JR ^ 

5o3e(2Z+i) 



(a;o,?o3)/ (a;i,X2,-||^p, 

d^^u {xi,X2, --^^p,t 
),Q(«o3) (a;o,^o3)/ (a;i,X2, 

d^^u {xi,X2, il^p, 
Recall that for any (P, Q) G N^, 

As,P,Q iXo,^o3)f (xi,X2, ^iP,tj =0 Vs > , 



dxidx2dtds 

(2.18) 



dxidx2dtds . 



(2.19) 



-4s,P,Q {Xo,S.o3)f {xi,X2,--^p,tj 

= A-2e(x„,eo3)/(a;i,a;2,-|jp,t) +^f+Ma;o,eo3)/(a:i,X2,~^p,i) Vs > , 

where 
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< ^ /r. S^^^l 1^ (C, r) I \a (0, 0, (4Q + 1) + + 2(,^hs, t + 2t/is, , 



Ar+M^o,Co3)/( 

< ^ 1^/ (e, r) I I a (0, 0, (4F + 3) - x,, - 2^,hs, t + 2Ths, s 



d^dr 



Now, for L > 0, Co3 e (2Z+1) and (a;o3,6) e hi f ] x Ko3 - 1, Cos + 1], we choose Q (^03) and 
-P (C03) large enough, for example 



o3) 



^(Co3)- P-^[(i + l) + 2(|eo3l + l)/^oi] 

in order that for any s S [0, L], /i € (0, ho] and (2:03, ^3) e [-|, f ] x [^03 - 1, ^^3 + 1], 



a (0, 0, (4g + 1) + xoz + 2^/15, 0, 
a (0, 0, - (4P + 3) ilfi p + a;o3 + 2^3/15, 0, s 



< 



< 



e ■4'' 



e . 



Consequently, 



?o3 6(2Z+l) 



{\dx3u{xi,X2,-p,t)\ + \dxsU{xi,X2,P,t)\) 



(f + 2Tha)^ 1 



and we conclude, using p.l8|l - p.l9|l . a classical trace theorem and l|2.8(l . that 



*^/o /fr,ur2)xR ^,P(C„3),Q(Co3) (2^0,^03) /(a;, t)5i.w(x,t)da;dids 

?o3e(2Z+i) 

Twf /r/r3 /|r|<A 



/ 1 



< e-ir L c^g idtu,0) {^Y ^/gJ^) 



d^dT 



dtds 



(t + 2Tha)^ 1 

e 4 (TTiyiiTT d^dTy/g {dtu,t) dtds 



d^dr 



(2.20) 



Now, we study to following boundary term 



ih 



^s,-P(«o3),Q(Co3) (2^0, C03) / (x, t) duu (x, t) dxdtds . 

^~XR , ^,0^ , 1^ 



eo3e(2Z+i) 



First, it holds 



As,P,Q {Xo,U) f {X,t) < J2 h" {Xo,U)f{x,t)\ 



(2.21) 



For ^o3 e (2Z + 1), we have 

A^{xo,U)f{x,t) 



< 



11 



Noticing 



we deduce that 

E |A^(x„,Co3)/(a:,i)| 

On another hand, we get 



d^dr 



(2.22) 



/txR ^ 

<ItI 



\t+2Ths\<\ 



(f + 2Th3)^ 1 

« T^^^i \d^u{x,t)\dxdt (2.23) 



^ /T/|t|<|2r/is|+V7rs+-U \dv,uix,t)\dxdt + c{l + hs)e bh y/g {dtu, 0) 



by cutting the integral over t S M into two parts and using a classical trace theorem. From (|2.21|) . 
(1223, it^^ and (EiHll, we conclude that 



JTxR X] As^p^Q{xo,^o3)f{x,t)d^u{x,t)dxdtds 
C„3e(2Z+i) 



< 



< 



JTxR X] X] l-4"(2;o,^o3)/(a;,i)l |5i'W(2;,^)l'^2;dtds 

C„3G(2Z+1) riGZ 

Co3e(2Z+i) 



/txR ^ 



(t + 2rfc»)^ 1 



\t\<2XhL + VhL + ^ 



< L 



i_ \d^u {x, t) I dxdt + c (1 + /is) e (9tM, 0) ) ds 

s 1/2 



(2.24) 



Finally, and I^TM imply 

*^/o^/anxR E ^,P(e„3),Q(«o3) K'^o3)/(a;,05!^"(2;,i)c;2;(^^c^s 

4o3e(2Z+i) 

1 /2 

< L c (A)^ (/^ /|t|<7(^)^ 1^-" '^^^O + ^ c (A)'' e-i y^lMlx/^ (9tu,0) . 

(2.25) 



2.9 The choice of A and L 

By (123, I2IIU, (|?^ and ^T^, we obtain, when / = d'^u, 

<c^y^T^V^(5tw,o) 

+c(^)^e-tg (zi,0) 

(2.26) 
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We choose A > 1 and L > 1 be such that ■j-^^= — y/h and -j-^-^ i^)'^ — in order that 

1/2 (2 2T) 

< cVT^x/^KO) (St?^, 0) + c (i)^ (/^ /|t|<7(i)^ 1^-^ '^^'^^) 7^^^- 

By replacing / by w and a by a sohition of (ids + ^ ~ ^t)) ^ (^^j t, s) = given by 

such that a{x — Xo, t, 0) = Uo [x, t/-\/2) ,we can argue in a similar fashion than above to show that 
T^r /oxE W2) (a;, t/\/2) u (x, t) u (x, i) dxdt 

1/2 (2 28) 

< cVh^/gKO)^/S [dtu, 0) + c (i)^ ( /|t|<^( 1 )^ \d.u (x, t)f dxdt) ^/W^) ■ 

Consequently, from p.2() . H2.27|l - H2.28() . we obtain that for any {xq} .^^ e uJZ, 
/axE Xa;' (a;) \a [x - x^, i, O) dtu (x, i) |^ dxdt 

< cy/h^/gj^^g {dtu, 0) + c (i)^ (/^ /|ti<7(i)'' 1^"" '^^'^O ^'^ V^T^ , 

and finally, by (|2.3|l . 

L„x(o,T) |9tw(x,i)|^dxdt 
< CoVh^G i% 0)VG (dtu, 0) + Co (i)'' ( /|(|<^(^)- |9,u(x,<)|'dxdt)'^V^K0) . 
for some Co > independent of m and h. 

3 Proof of the main result 

The choice Wo = {—mi + ro, mi — Tq) x (— m2 + To, m2 — To) x (— f , f ) allows to invoke the geometrical 
control condition. Indeed, observe that any generalized ray of — A parametrized by i € [0,T], for 
some T < +oo meets Wo U w. As a result, we have the following observability estimate. There exists 
C > such that for any initial data {uq, ui) E Hq (f2) x (f]), the solution u of (|1.2|l satisfies 

II(wo,wi)|Ih1(j7)xL2(o) < C ^ \dtu{x,t)f dxdt + J \dtu{x,t)f dxdtj . (3.1) 

From Theorem 2, the solution u of 11.2|l satisfies the following interpolation inequality. There exist 
C > and 7 > 1 such that for any h G (0, ho], 

1/2 

LJo \^tu{x,t)\'^dxdt <C{j^y Q^J^^^^^^iy\d^u{x,t)fdxdtj IIK,Ul)||jji(o)xL2(i2) (3 2) 

+CVh ||(M0,Ml)||ff2nHi(o)x//i(O) ll("0,'«l)|l_f/i(n)xL2(o) • 

It is now known | Li l using the multipliers technique that the normal derivative of the solution of the 
wave equation with Dirichlet boundary condition satisfies the following inequalities. Let ip G (0) 
be such that i/) = 1 on T, then there is c > such that for any T > 0, 

||9i,u||^2(xx(-T,T)) - II^^IIh1(wx(-2T,2T)) 



' (ll^t'"llL2(wx(-3T,3T)) + I1'"I1l2(wx(-3T,3T))) 
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Consequently, from l|3.1|) - (|3.2|) - (|3.3|) . a translation in time, we obtain that the solution u of (|1.2|l 
satisfies 



+CVh \\{uo 



As dtw, i.e., the derivative in time of w solution of Hl.ll) . can be seen as a solution of the wave with 
a second member —a {x) dtw, (|3.4(l implies with a usual decomposition method, knowing a > on to 
and £ {w, 0) < cS {dtw, 0), that there exist C" > and 7 > 1 such that 

£ (w, 0)+£ (dtW, 0) < C {j^y Jo^^^'' (a (x) \dfw {x, t)| V a (x) \dtw {x, dxdt 
+C'h £ {dfw, 0) Vh G (0, ho] . 

This later estimate is clearly also true for any h > ho and some constant C" > 0, thus we can choose 

1 £{w,0)+£{dtw,0) 
~ (C + C") £{d!w,0) 

in order that there exists C > such that 



^ \£{w,0)+£(dtw,0) J JoJc 

By a translation in time, we obtain that 



e 



£iw,s)+£{dtw,s) ^ r< f £{d^wfi) Y I- £(™,s)+£Ot»,a) I a(x)\a^w{x,t)\^+aix)\dtw(x,t)\'' ^ 

£id'iw,0) - ^ {£{w,s)+£{dtw,s) J JnJs £{dfw,0) "'^"'^ 

^ r ( ^jSt-^.o) V / £{w,s)+£{dtw,s) \ ( £{w,t)+£ (dtw, t) \ 

\£(w,s)+£(dtw,s) ] \ £(d'iwfi) i I £(d'twfl) J , £(a'i^,o) 



Applying Lemma in Appendix B to 

£ {w, s) + £ {dtW, s) 



£{dfw,0) 



where ct > is taken in order that !F is bounded by one, we get that there are C > and S > 0, such 
that for any t > and any initial data ({uq, wi) € (Q) n Hq (Q) x (fl) n Hq (fl), the solution w 
of 

dfw — Aw + a {x) dtW — in x M 

w = Aw ^0 on afi X M (3.5) 

(w(-,0) ,5tw(-,0)) ^iwo,wi) in , 

satisfies 

^ |V(9tw(a;,i)|^) da; < ^ ||(wo>i)||ff3(n)xij2(o) . (3.6) 

Now, it is known ' Li how to deduce from H3.6|l that there are C > and S > 0, such that for any 
t > and any initial data (wq, wi) € i?^ (il) n Hq (ft) x Hq (O), the solution w of Hl.l|l satisfies, 

\dtw{x,t)f + \Vw{x,t)f^ - F "(^0'"'i)"^f"(S2)xHMf^) • 
Indeed, let wq G i?^ (fi) n Hq (il) be such that Awq — wi + awo E Hq (fl). In particular, it comes 

||wo||^3(o) < cllwi +awo||^i(f^) . 
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Now, notice that w solves 

dtw {x, t) — wi {x) — /p Aw {x, t) dr + a {x) {w {x, t) — wq (x)) = , {x, t) & Q x , 

w = on 9f] X R+ , 
w (•, 0) = wq, 9(W (•, 0) = lui in , 

and as a result, w {x, t) — w {x, r) dr + wq {x) solves (|3.5|) with dtW (-,0) — wq and satisfies H3.6|) . 
And we conclude that 

(\dtw{x,t)f + \Vwix,t)f^ dx <^ (c||wi +awo||^i(s:2) + ||wo||^2(o)) • 
That completes the proof. 



4 Appendix A 



The goal of this Appendix is to prove, with the notations of the above sections, the two following 
inequalities, 



/oxR (2^)^ /r3 i|r|>A e<'=^+^^'^ (pu (C, t) d^dr u (x, t) dxdt 

+ /nxR"o(^'*) (2^/R3/|r|>Ae''*''^^*''^ VdtU{i,T) d^dTu{x,t) dxdt 



and 



'|r|<A 

Proof of (|AH) . Introduce 

R{9) = 



|r|<A 



d(,dT <c\^\g (u,0) 



/nxR«°(2^'^) ^*((2i)^/R3/|r|>Ai <^(e,T) dCdr) u{x,t)dxdt 

~ In.Mdt{aou{x,t)) {^^ J^^^^.i^ e^'- [J^e-^<>-^g{x,9) dO] dr) dxdt 

It follows using Cauchy-Schwartz inequality and Parseval identity that 



R{g) < \dt {aou (x, t))\ (2^y^/|,|>, ^dry^/j, 1/^ e-^»-^g (x, 0) d0| 



^ /nxRl^t(«o"(2^'^))l (^\//|r|>A7?'^^\/27r4|(^.g(x,6l)|'^ 

^ /nxR l^t (2^' ll'^ff (2^' •)IIl2(r)) rfa;dt 

/Rll^«(ao")(-'*)llL2(a)C;^ll¥'5(a;,-)llL2(iixR) • 
Since we have the following estimates 



dO dxdt 



ll¥'"(a;,-)llL2(nxR) + lk9t^"(2;r)|L2(oxR) ^ \// Rg ^''/n |M(a:^,t)|^dxdt 

+ \//Re~^*'/o |9?u(x,t)|^dxdi 
< Vvl^ (c^^TM) + \/e (5tM, 0)) 



(Al) 



(A2) 



/jg ||5t (aow) (■,t)||^2(f2) di < /r -\/ l^tOoU (x,t)|^ dxdt + Js^\J \aodtu (x, dxdt 

< /r Jt£^ vtK^^'^^* + -^K \[kAd^u(^^^dxdt 
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we conclude that 

R{u) + R (d^u) < c^^/g {u,o)^g{dtu,o) . 

That completes the proof of HA1(I . 

Proof of IIA2|l . We estimate /r3 d£,dT wheie g solves df g — Ag ^ in il x M.. Recall 
that f{-,t) — xo-o ('ji) G C^' (f2). So, using Cauchy-Schwartz inequality and Parseval identity, 



/m3 J^r|<Al'^(^'^)l ^^^^ =lR^I\r 



i±iii 



^\'P9iC,r)\ d^dr 



d^dr 



d^dr 



<^'A//r3/m<a (l-A)(</'5)(C,r) 



d^dT 



Observe that 



A {<fg) = ifAg + 2V(/?Vg + A(pg 

= ifid^g + 2ip {xo, t) Vip (x, 0) Vg + Aipg 

= 92 (^^) _ 2dt (dtipg) + {Aip + dfif) g + 2ip (xo, t) Vip {x, 0) V.g 



As a result, 



|r|<A 



fc=0,l,2 



d^dr 



+cVA^/r3 /|,|<;, ^ {xo, t) (x, 0) Vg r) 
In particular, when g = u we obtain, using Parseval identity 

+cVA 11^ (xo, t) Vip {x, 0) Vu|| i2(oxRi 



< c 



VAA^+cVAi) /ije-^*' (^|u(x,t)|'dx) dt + cVX^^/G (u, 0) 



On another hand, when g = dfu and using the fact that |r| < A, 



/|r|<A 



\ /r3 /|t|<A ^ 

y fe=o, 1,2,3 




2 

d^dr 




/r3 /|t|<A S 
fc=0,l 




2 


Iw^ I\t\<X S 
J fc=0,l,2 


V (^o,^*)^^' (x, 0) Vu (e, r) 


2 

d^dr 



u,0) . 

We conclude that there exist c > and 7 > 1 such that for any h e (0, ho] and A > 1, 



|^(^,r)| d^dr 

/|r|<A 

That completes the proof of l|A2|l . 



R3 J\t\<\ 
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5 Appendix B 



Lemma.- Let J- be a continuous positive decreasing real function on [0, +00) and bounded by one. 
Suppose that there are four constants ci > 1 and C2,/?, 7 > such that 

Then there exist C > and 6 > such that for any t > 0, 



Proof .- Let t > 0. If (^)^ < 7 then J^(s) < If } < (^)^ then (^)'' + s<t + s, 

thus {J^{t + s))<J' ((f^) ^ + ^^'^ therefore 

jr(s)<(ci)3TT{jr(s)-(jr(f + s))}wT . 

Consequently, 

^(s) < (ci)^ {J^{s) - (^(t + + -0- Vs,i > . 

Let 



^Cii^ (3 + 1 _j_ C2 



1 

If .F(s) - + < j^, thcnTis) < (f^)^ + 70t and thus ibt {t + s) {t + s) < 1. If 
^ < J- (s) - ( J- + s)), then ^ ( J- (s)) < ^ < J- (s) - + s)). Therefore 



iPt{t + s)J^{t + s) <^jr(s)^t(t + s) = Vt(s)J^(s) ^^^27r 
<Vt(s):^(s) 

by using the decreasing property of C. ' — > ^"^''^ ■ We have proved that for any s,t > 0, we have either 
'ipt {t + s) {t + s) < 1, or ipt {t + s) {t + s) < tl^t (s) (s) . In particular, we deduce that for any 
i > and n e N\{0} , either 

^Pt{{n+l)t)Tiin+l)t) < 1 
or Vt {{n + 1) ((n + 1) < V't (rii) ^ ("0 • 

Then inductively, it implies that 

ipt {{n + 1) i) ((n + 1) < max (1, Vt (t) (t)) = 1 • 

Hence for alH > and n e N \{0} , 

1 

Ci \ f' + i C2 



We choose n such that n + 1 <t <n + 2 and we obtain that for alH > 2, 



The desired result now follows immediately. 
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